Response to the Comment on "Excitons in Molecular Aggregates with Levy Disorder: 
Anomalous Localization and Exchange Broadening of Optical Spectra" 
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In our Letter [l|, we predicted exchange broadening 
and a blue shift of the absorption spectrum with increas- 
ing disorder strength a for heavy-tailed disorder with a 
stability index a < 1. We explained these findings by de- 
riving 'conventional' scaling laws for the width of the ab- 
sorption spectrum and the delocalization length. In addi- 
tion we predicted new features in the density of states and 
the absorption spectrum, as well as non-universality of 
the localization length distribution, which we attributed 
to the appearance of 'outliers' in the disorder, i.e. sites 
with energies outside of the bare exciton band. 

The main points of the Comment are that the out- 
liers introduced by heavy tails in the disorder distribution 
(i) do not lead to deviations from the conventional scal- 
ing law for the half width at half maximum (HWHM) 
of the absorption spectrum and (ii) do not lead to non- 
universality of the distribution of localization lengths. 
We show below that the findings reported in our Let- 
ter [H are correct and that the wrong conclusions of the 
Comment Q arise from focusing on small a values. 

The conventional scaling law for the HWHM 
(= FWHM/2) in the general case of symmetric a-stable 
Levy distributions reads 
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HWHM - J(cr/J)2"/(i+"). 
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Its derivation was given in our Letter [ij already, gener- 
alizing from the special cases of Gaussian (a = 2) and 
Lorcntzian (a = 1) disorder |3l to arbitrary values 
of a. In the Comment Q], this derivation is repeated 
without any additional insights (last paragraph of the 
Comment In the derivation, segmentation due to 

outliers (defined by us as sites with deviations from the 
average energy larger than 2 J, acting as barriers) as well 
as finite size effects are not taken into account. For a 
values of interest for many molecular systems (cr ^ J), 
deviations from the conventional HWHM scaling due to 
outliers have not been seen for a = 1 or 2. However, for 
a « 0.5 and smaller clear deviations are visible in this a 
regime. This is due to the fact that for the same cr, out- 
liers become more abundant for heavier tails (smaller a). 
By equating N* and iVgog from Eqs. (3) and (4) of our 
Letter one can estimate the value of a above which 
segmentation may be expected to play a significant role 
(indicated as the shaded region in Fig. la). 

To demonstrate the deviations from the conventional 
scaling at smaller a values, we plot in Fig. 1 the HWHM 
as a function of a for a ~ 0.5 (b) and a ~ 0.4 (c). The red 
lines are power-law fits imposing the conventional scal- 
ing law. The case a = 0.5, panel (b), was also given in 
our Letter [l[ but now has many more data points (sym- 
bols), also compared to the Comment [ij. For a < 0.4 J 
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FIG. 1: (a): "phase diagram" showing for which (a, cr) com- 
binations segmentation is expected to become visible. (b,c): 
HWHM of the absorption peak as a function of a for (b) 
a = 0.5 and (c) ot = 0.4. Data points obtained from numer- 
ical simulations (A*" = 200, 10^ realizations), red curve rep- 
resenting a fit corresponding to the conventional localization 
mechanism (note the deviations from the fit). 



the conventional power law fit works well, agreeing with 
the observation in Ref. [H and repeated in the Com- 
ment; for larger cr, the conventional scaling clearly breaks 
down. Moreover, the data points show kinks at which the 
HWHM suddenly jumps. This is not noise, but may be 
unambiguously attributed to structure in the high-energy 
wing of the absorption spectrum which grows in for in- 
creasing cr as a consequence of segmentation While 
both effects are visible already at a = 0.5, they become 
more pronounced for smaller a (panel c). This new data 
also shows that these features cannot be adequately rep- 
resented by a power law. Finally, we mention that the 
straight line in Fig. 2 of our Letter [l| indeed is meaning- 
less without the information how the monomer HWHM 
scales with cr; this has no effect on the conclusion of that 
figure, namely the discovery of exchange broadening. 

Next we demonstrate the non-universality of the local- 
ization length distribution due to the segmentation by 
outliers. While the data used for this purpose in our 
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FIG. 2: Distribution of normalized localization lengths for 
small disorder values (a) and intermediate disorder values (b). 
In (b) the curves have different vertical offsets to better dis- 
tinguish their differences. Notably, additional structure ap- 
pears at the low localization length side, corresponding to 
segmentation-limited wave functions. N = 200, a = 0.5. 
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Letter [l[ (Fig. 3 and the a-dependence of SNioc/Nioc) 
were correct, they were not a good choice, especially not 
because they were discussed in the context of Ref. [1], 
where the energy interval in which states are considered 
was scaled with the width of the absorption band. In 
Fig. 2, we show that also for an energy interval which 
is "properly" scaled, the non-universality becomes visi- 
ble. This figure displays for a = 0.5 the distribution of 
normalized (to the average) localization lengths, obtained 
using the same scaling of the width of the energy interval 
as in the Comment Q (cf. Ref. Q)- In Fig. 2(a) small a 
values arc considered (cr < O.IJ, as in the Comment 0]). 
As expected from Fig. 1(a) segmentation is not relevant 
then and the normalized distributions are nearly identi- 
cal for all (J < O.IJ. Yet, already for cr ~ O.IJ deviations 
become clearly visible in the appearance of several small 
sharp peaks at small localization length. Upon increas- 
ing cr (Fig. 2(b)), these deviations become more and more 
pronounced and the non-universal character of the dis- 
tribution is clear By considering the wavefunctions, 
the origin of these peaks can be traced back to the seg- 
mentation mechanism. We note that also for the fixed 
energy interval used in our Letter [l[ , these signatures of 
segmentation and non-universality can be found. 

In conclusion, we have shown that the claims put forth 
in the Comment are not generally valid, but are instead 
an artifact of focusing on a regime where segmentation 
effects are expected to be small. In the Comment it is 
stated that conventional scaling may break down at high 



cr values, only because Nioc approaches unity. This is a 
trivial effect and does not explain the additional structure 
observed in the localization length distributions. As we 
have demonstrated, the breakdown found by us is more 
subtle and really is intimately related to outliers. 
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